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Abstract 

Optimal observables are known to lead to minimal statistical errors on param- 
eters for a given normalised event distribution of a physics reaction. Thereby 
all statistical correlations are taken into account. Therefore, on the one hand 
they are a useful tool to extract values on a set of parameters from measured 
data. On the other hand one can calculate the minimal constraints on these 
parameters achievable by any data-analysis method for the specific reaction. 
In case the final states can be reconstructed without ambiguities optimal ob- 
servables have a particularly simple form. We give explicit formulae for the 
optimal observables for generic reactions in case of ambiguities in the recon- 
struction of the final state and for general parameterisation of the final-state 
phase space. 
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1 Introduction 



At a future e+e~ linear collider (LC) like TESLA [1] or CLIC [2], apart from the 
possible discovery of new particles, electroweak precision measurements of various 
observables will be an important task. In this way effects due to new physics at 
a scale far beyond the produced cm. energy may be detected. The e+e", the 77, 
the e'~7, the e~e'~ and the Giga-Z mode can reveal complementary aspects. For in- 
stance, consider the thrce-gauge-boson vertices jWW and ZWW, which are highly 
restricted in the Standard Model (SM), see for instance [3]. In a form-factor approach 
with the most general Lorentz-invariant parameterisation these vertices are described 
by 28 real parameters if one allows also for imaginary parts [4]. In many extensions 
of the SM one obtains deviations of these triple-gauge-boson couphngs (TGCs) from 
their SM values, for an overview of the literature see e.g. Sect. 1 of [5]. At present 
the experimental constraints on these couplings from LEP are rather weak and many 
couplings have yet to be measured [6-8]. At a future LC all TGCs can be measured 
with much higher precision, see [9-11] and references therein. Both longitudinal and 
transverse beam polarisation are considered to be feasible [12] and will be advanta- 
geous in many cases, see for example [13]. For most TGCs longitudinal polarisation 
is found to be the right choice [5], whereas one coupling can only be determined with 
transverse polarisation [14]. In a gauge invariant effcctivc-Lagrangian approach for 
the gauge-boson sector new-physics effects can be parameterised by ten anomalous 
couplings [15,16]. In such a framework, after spontaneous symmetry breaking also 
gauge-boson-fermion interactions and gauge-boson masses are modified by the anoma- 
lous couplings. Therefore constraints can be derived from observables measured at 
LEPl, SLC and LEP2. Also in this case the couplings can be determined with much 
higher precision at a future LC [16]. 

In all cases it is important to know how sensitive the event distribution of a 
certain reaction is to a given set of parameters, independent of the method that will 
be used in the future to analyse the data. Moreover statistical correlations should be 
taken into account, unless they are very small, in order to provide realistic results. 
Frequently one is dealing with a situation where the normalised event distribution 
is — at least approximately — a linear function of the parameters to be measured. The 
method of optimal observables [9, 17] is most directly applied for the case of a strictly 
linear function. In most parts of this paper we assume that the event distributions 
are to good approximation linear functions of the parameters to be measured. By 
means of optimal observables one can then compute the maximum constraints on 
this set of parameters for a given reaction and given event number while taking into 
account all statistical correlations. Apart from being a useful tool for theorists to 
estimate the sensitivity of a reaction to a set of parameters optimal observables have 
been used in experimental analyses to extract the parameters, for instance TGCs [6] or 
CP violating parameters in e^e" — > r+r^ [18]. But also methods to handle non-linear 
functions of the parameters to be measured have been devised [10] and successfully 
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applied in experimental analyses, see e.g. [8,19]. 

In case of many parameters with correlated errors, the sensitivity to different di- 
rections in parameter space is often not easy to survey. Here optimal observables 
have the advantage that discrete-symmetry properties of the differential cross section 
can be exploited in order to eliminate correlations between couplings. For example 
the TGCs are classified into four symmetry groups, and couplings from different sym- 
metry groups can be measured without correlations when optimal observables are 
used [9]. However within each symmetry group couplings are in general correlated. 
In [10] it is shown that the optimal observables are unique up to linear transforma- 
tions. 

Here we construct explicitly the optimal observables for generic ambiguities in the 
reconstruction of the final state. Such ambiguities can be either continuous, that is 
the number of measured variables is smaller than the number of variables required 
to specify the final state, or discrete, that is a reconstructed event can originate from 
two or more final states. 

This work is organised as follows: In Sect. |21we resume the properties of optimal 
observables as given in [9]. In Sect. El the optimal observables for generic ambiguities 
in the reconstruction of the final state are discussed. In Sect. |3]we recall that the 
method can be applied iteratively to determine parameters on which the differential 
cross section depends non-linearly. We present our conclusions in Sect. El 

2 Optimal observables 

In this section we give a resume of the definition and properties of optimal observ- 
ables. As it is convenient to have an illustrative example in mind we shall discuss 
the problem to measure anomalous contributions to the differential cross sections for 
the reactions e^e~ — > WW and 77 — > WW. But we emphasise that our considera- 
tions are neither restricted to anomalous couplings nor to particular reactions. The 
optimal-observables method can be applied to any reaction where the differential cross 
section depends on a certain number of small parameters, which we generically denote 
by hi and which are to be estimated. Furthermore, the method can be generalised 
to the case where the parameters to extract are not necessarily small [10]. Consider 
now the reaction e+e~ WW or 77 WW and let us assume that we describe it 
in the framework of the SM with the addition of small real anomalous constants hi. 
In an experiment one measures the differential cross section 

S{4>) = da/d^, (1) 

where (p denotes the set of all measured phase-space variables. For instance, the 
fully differential cross section of the process e~^e~ WW without transverse beam 
polarisation depends on five angles, which are in this case specified by 0, see [5]. In 
the same way the spin-averaged fully differential cross section of 77 — > WW with 
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fixed photon energies depends on five angles, see [20]. In case of e^e~ — > WW as 
treated within a form-factor approach in [5] the anomalous parameters hi are the 
28 anomalous TGCs that parameterise deviations at the 'fWW and ZWW vertices 
from the SM. In case of 77 WW as treated within an effective-Lagrangian approach 
in [20] the hi are the ten coefficients of certain dimension-six operators. Here non- 
zero anomalous couplings do not only lead to anomalous three- and four-gauge-boson 
couplings but also to deviations of the gauge-boson-fermion couplings and of the 
gauge-boson masses from their SM values. We distinguish between the information 
from the total cross section a = J da and from the normalised distribution S/a of 
the events. Here we only investigate how well anomalous couplings can be extracted 
from the latter. It is possible to obtain constraints on these parameters also from the 
measurement of a, see Sect. 3.1 of [10]. Those considerations remain unaffected in 
the presence of ambiguities in the reconstruction of the final state. 
Expanding 5* in the anomalous couplings one can write 

S{<P) = So(0) + J2 Sim + Oih^), (2) 

i 

where So{(j)) is the cross section in the SM and the Sii{(f)) give the first-order modifi- 
cations due to the anomalous couplings. We assume So{(j)) and Sui^cf)) to be calculated 
from theory. Note that the variables need not specify the final state completely. In 
the analyses [5, 14] of e'^e~ WW with one W decaying leptonically and the other 
one into two hadronic jets it was assumed that the jet charges cannot be identified, 
which results in a two-fold ambiguity. In such cases S{(j)) is not the fully differential 
cross section — in our example it is the sum over two final states. In a case with 
ambiguities of this or more involved kind it is often not straightforward to calculate 
5*0(0) and Sii{(j)). But these quantities must be known explicitly in order to construct 
the optimal observables. This problem is our main concern in this paper and will be 
addressed in the following section. 

We now give a short resume of the optimal-observables method. One way to 
extract the parameters hi from the measured distribution Q is to look for a suitable 
set of observables Oi(0) whose expectation values 

E[Oi] = - U<PS{<P)0,{<P) (3) 
a J 

are sensitive to the dependence of 5'(0) on the couplings hi. To first order in the 
anomalous couplings we have 

E[Oi\ = E,[Oi\ + c^Jh, + 0(/i2), (4) 
3 

with 

Eo[0^] = - fd<pSoi<P)0,{<P) , (5) 
an J 
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Cij 



^ 0,{<P) fd<p So{<P) 0,{<j)) , (6) 



(To J CTq 

ao = /d0So(0), (7) 



di, = yd0Si,(0). (8) 

Here -Eo[C^j] is the expectation value for zero anomalous couplings, and Cij gives the 
sensitivity of to hj. Solving for the set of the hj we get estimators for the 

anomalous couplings, whose covariance matrix is given by 

V{h) = j^c-'V{0){c-y, (9) 

where we use matrix notation. Here is the number of events, and 

V{0),, = - I d05o(0) - E^Pi] EoPj] + 0{h) (10) 

(Jo J 

is the covariance matrix of the observables, which we have expanded around its value 
in the SM. As observables we choose 

om = (u) 

From (jH)) and (fTUI) one obtains for this specific choice of observables 

ViO) = c + Oih), (12) 

and therefore 

Vih) = ^c-' + 0{h). (13) 

From (|12j) we see that for the observables (|lip c is a symmetric matrix because V{0) 
is symmetric. The observables (jlip are "optimal" in the sense that for hi ^ the 
errors ()13p on the couplings are as small as they can be for a given probability dis- 
tribution, see [9]. For details on this so-called Rao- Cramer- Frechet bound, see for 
example [21]. Apart from being useful for actual experimental analyses, the observ- 
ables (jllj) thus provide insight into the sensitivity that is at best attainable by any 
method, given a certain process and specified experimental conditions. In case of 
one parameter this type of observable was first proposed in [17], the generalisation 
to several parameters was made in [9]. Moreover, it has been shown that optimal 
observables are unique up to a linear reparameterisation [10]. We further note that 
phase-space cuts, as well as detector efficiency and acceptance have no influence on 
the observables being "optimal" in the above sense, since their effects drop out in the 
ratio (fTT|) . This is not the case for detector-resolution effects, but the observables (fTT|) 
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are still close to optimal if such effects do not significantly distort the differential dis- 
tributions Sii and 5*0 (or tend to cancel in their ratio). To the extent that they are 
taken into account in the data analysis, none of these experimental effects will bias 
the estimators. 

When a set of events is calculated with a Monte-Carlo generator and the events 
are then reconstructed like in an experimental analysis the optimal observables (|TT|) 
are obtained directly in the measured variables 0. However for theoretical studies, i.e. 
to estimate the sensitivity of a certain reaction to anomalous couplings, and also for 
experimental analyses often the analytical expressions of the optimal observables (fTT|) 
are required. 

Frequently there are ambiguities in phase space, that is to one value of the mea- 
sured kinematic variables (p there correspond two or more distinct final states (dis- 
crete ambiguities) or a bunch of final states (continuous ambiguities). The calculation 
of So{(j)), Sii{(j)) in (0) and in particular of the observables Ci(0) in (fTT|) has then to 
be done with some care as will be shown below. 

3 Phase-space ambiguities 

In principle there are plenty of possibilities to parameterise a final state in a reaction 
uniquely, for instance the usage of angles or Cartesian coordinates, different choices of 
reference frames etc. In an experiment one may either be able to specify a final state 
of an event uniquely or only with certain ambiguities. Here we discuss in detail the 
case of discrete ambiguities, that is for each event one only knows that it belongs to 
a group of two, three or more final states. We also mention how to handle continuous 
ambiguities, i.e. the case in which the number of measured variables is smaller than 
the number of variables required to specify the final state. An example of a discrete 
ambiguity is the two-fold one of the semileptonic final states in e^e~ WW or 
in 77 — > WW with fixed cm. energy of the two-photon system. Here one usually 
assumes that the two hadronic jets cannot be associated unambiguously to the quark 
and antiquark. Another more involved one occurs in the reaction 77 — ^ WW when the 
photons each obey a Compton spectrum. Here, in addition to the ambiguity above, 
another two-fold one arises from the reconstruction of the neutrino momentum. This 
case is considered in [20]. The optimal observables (jllll are now to be calculated in 
the presence of such ambiguities. 

We start from a particular set of phase-space variables x fhat specify the final 
state uniquely. The differential cross section in terms of these variables we denote by 



The cross section for another choice of variables 0, that may lead to the above ambi- 
guities, is then given by 



T{x) = do-/dx. 



(14) 




(15) 
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The function F expressing the relation of to x niay take the same value for different 
values of x, that is for a given the equation 



Fix) = <P (16) 

may have several solutions Xk = Xk{4') with k = 1,2, . . .. In general, the number of 
solutions to (fTHIl may vary with (p. If (p are the coordinates that can be measured 
of an event x, the set of final states Xk consists of x itself as well as all final states 
that cannot be distinguished from x by a measurement of (p. Notice that in ()15|) we 
have assumed that the number of components of the vectors and x are the same. 
In general there can be also continuous ambiguities, that is the fully differential cross 
section is specified by a larger number of variables than those that can actually be 
measured. However this does not lead to any further complicacies in the context of 
optimal observables. In fact, if the fully differential cross section is T(x, where the 
final-state variables ^ cannot be measured, the generalisation of (fT3j) is 



5(0) = / dxde S{F{x) - 0) T(x, 0- (17) 

If we define 



T(x) = / d^T(x,0 



we again obtain ()15|1 . We can thus apply all formulae in the remainder of this section 
also in case of discrete plus continuous ambiguities. We remark that our analysis also 
works if one or more of the phase-space variables take discrete values as is the case 
for instance for spin indices. For these variables integrals have to be substituted by 
sums and ^-distributions by Kronecker symbols. 

An integration and summation over part of the phase-space variables as in ()17j] is, 
of course, performed when one considers inclusive cross sections. Thus our discussion 
covers also this case. Clearly, then the normalisation integral / d0 S{(j)) gives the 
cross section times the corresponding multiplicity and in our formulae a has to be 
read in this way. 

Coming back to p5j] we have 

S{<l>) = J2\Jk\-'T{xk{<P)) (19) 

k 

where 

Jfc = det|^(Xfc(0)) (20) 
OX 

is the Jacobian determinant taken at point Xk- If F is invertible, there is only one 
term in the sum for all and (fT^ simplifies to 



5(0) 



t(f~\<P)). (21) 
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We expand the differential cross section: 

T{x) = To{x) + Y.Tu{x)h. + 0{h'). (22) 

i 

It follows 

Si(l)) = Soi(l)) + J2Simhi + 0{h^), (23) 

i 

where 

SM = T.\M"'n{xkm. (24) 

k 

Su{(t^) = Y.\M''Tu{Xkm- (25) 

k 

Note, again, that the number of terms in the sums and can vary with 0. If 
(p — but not necessarily x — are coordinates that can be measured we have to define 
the optimal observables from the expansion of 5(0) in (|^^: 

Om = (26) 

In the specific case where F is invertible going from x to is a mere change of 
coordinates and we obtain the same optimal observables using either set of variables: 



Toix) 



(27) 



If there are ambiguities in the reconstruction but if we have the same Jacobian J = Jj. 
for all k (which may nevertheless depend on 0), J cancels in the numerator and 
denominator of the observables (|26p: 

Oi{(f)) = — -— . 28 
EkToiXkW) 

This is the case e.g. for the reaction e^e^ — > WW, where one W boson decays into 
a quark-antiquark pair and the other one into a lepton pair, if the charges of the two 
jets in the final state cannot be identified, see [9]. If this is not the case we must use 
the general expressions (j2H) to pHjl . 

The covariance matrix of the observables (1261) is now 



Vm, = - fd<pSo{<P)0,{<P)0,{<P)-^^ + Oih), (29) 
an J A an 



ao J A ctq 

where 



ao = / d0 5o(0)= / dxTo(x), (30) 

J A Jb 

an = / d0 S'i,(0) = / dx Tii(x), (31) 



B 
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F 



Figure 1: Definition of integration areas. 



and the full kinematically allowed integration regions in the coordinates and x ^-re 
denoted by A and B, respectively. The integrals (Tq and an can be performed in 
either coordinates. Using x iio knowledge about ambiguities in the reconstruction is 
necessary. The first term in the expression of V{0) needs special care. We divide the 
integration region A into parts An with n = 1,2,..., such that for (p & An there are 
n solutions Xk to (fTBj) . The domains of B corresponding to the An we denote by Bn, 
see Fig. [TJ We subdivide Bn into n appropriate regions Bnk with k = 1,2, . . . ,n, such 
that Xk £ Bnk- This subdivision is certainly not unique. We have 




(32) 





Sii{(p)Sij{(f>) 
5o(0) 

Sii{(p)Sij{(p) 
'^^ ^7T\ 



Expressed in terms of integrals 



over X we get 




+ E /„ dx|J(x)l 



n>2^"Pn 



Su{F{x))S^,{F{x)) 
So (Fix)) 



(33) 



and also 



H^, = f dx 



Toix) 



n>2 



Su{F{x))S^,{F{x)) 
So (Fix)) 



(34) 



where we have from (j21I), ()25|) for x G -Bj 



n 




(35) 



k 
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Su{F{x)) = j:\j{xk{F{x)))l T,,{xk{F{x))), (36) 

k 

J{X) = det^ix). (37) 

One of the values XkiF{x)) in fjH3j) . (jHUj) is, of course, identical to x- dSSl) one can 
choose for each n any natural number p„ with 1 < p„ < n. These choices correspond 
to different parameterisations of the integration regions An but leave the integrals 
unchanged. Therefore one may sum over all possible choices and divide by n, which 
leads to ()34p . The quantities Hy may be calculated either in the form (jHHj) or ()34|) . 
Notice that the form has the advantage that one only has to know where in the 
integration region for x there are how many solutions to (fTB|) . but one does not have 
to specify Bni, Bn2, etc. In certain cases the integrals for n > 2 in or (jH^ may 
be simplified. For example let A'^ with > 2 be the part of An where the Jacobians 
J{Xk{4')) are the same for all k. The Jacobian in this region may nevertheless depend 
on 0. The region of An where they are not the same for all k we call A'^. The 
corresponding regions of Bn are denoted by B'^ and B'^. We write the integrals 
in (jnH) as 

/ dx= f dx+ f dx. (38) 
Jb„ Jb'„ Jbh 

Then, in the integrals over 5^ the Jacobian cancels and we obtain the following 
expression for the integral in the covariance matrix (j29|) : 

rr _ f . Tu{x)Tlj{x) .^QS 

""'^ - Ux) ^''^ 

_1 r Y.lTu{xk{F{x)))Y:iT,,{xi{F{x)) 



^2"^^^'^ Ein(xm{F{x))) 

^2 ^ 



So (Fix)) 



with So{F{x)) and Sii{F{x)) as in and (IHUj). respectively. 



4 Iterative analysis in the non-linear case 

In this section we recall briefly that the use of optimal observables is not restricted to 
a phase-space distribution depending only linearly on parameters hi which are to be 
estimated. In other words, the higher-order terms in hi in (j21) can be handled. This 
has been discussed extensively in [10]. We recall here only one practical procedure 
one can follow. 
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Suppose we have a theoretical expression for the differential cross section Q which 
can be expanded in the hf. 



S{<j)) = SM + Y Su{<P) h^ + T. •52^,(0) h,h, + .... (40) 

i ij 

With given data an estimate of the hi has to be made. The procedure proposed in [10] 
is then as follows. In the first step the terms of second and higher order are neglected 
and one follows the procedure of estimating the hi by the optimal observables 
Suppose that this gives as estimates for the parameters the values hi. In the second 
step one sets 

hi = k + K (41) 
and substitutes this for S{(f)) in (jlUj) . This gives 

S{<P) = So{<P) + Y: Su{<P) K + Y S2i,{<P) Kh', + . . . , (42) 

i ij 

where 

So{^) = So{^)+YSl^{<P)h^ + Y.S2^,{^)h^hJ + ■■■^ (43) 

i ij 

Sli (0) = Sui<P) + 2Y S2r, (0) /i, + . . . , (44) 

i 

and so on. 

Now one applies the optimal-observables method to estimate the h[, neglecting 
terms of second and higher order in the h'i in fl42j) . The new optimal observables are 

Om = (45) 

Let h'i be the estimates for these parameters obtained in this way. The improved 
estimate for the original parameters is then hi + h[. 

This procedure can be iterated. It was tested in [19] in a Monte-Carlo study for 
the analysis of TGCs at LEP2. Parameters hi — not necessarily small — were assumed 
and Monte-Carlo data generated according to the corresponding distribution (jin|) 
which for this case contained linear and quadratic terms in the hi. The non-linear 
optimal-observables analysis was then performed as mentioned above. It turned out 
that after a few, typically three, iterations the input values for the hi were obtained 
back within their correct statistical errors. 

Clearly, the phase-space ambiguities are to be treated in exactly the way discussed 
in Sect. EJalso for the case of a non-linear analysis. 
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5 Conclusions 



For electroweak precision measurements at a future LC it will be important to check 
the validity of the SM (or perhaps another theory) with the highest possible preci- 
sion. To this end optimal observables are a convenient means because parameters 
can be determined with minimal errors as allowed by a theorem from mathematical 
statistics without neglecting correlations between any of them. Such observables have 
for instance been applied to the reaction e^e~ WW in [9, 10] with a form-factor 
approach to the 'yWW and ZWW vertices. Effects of beam polarisation to the same 
process were analysed in [5,14]. However in the mentioned studies the final state was 
assumed to be known either exactly or up to a two-fold ambiguity of a very simple 
type. It was possible to add the respective terms of the differential cross section 
in order to construct the optimal observables. In this paper we have discussed the 
calculation of optimal observables for the case of an arbitrary reaction with generic 
ambiguities in the reconstruction of the final states. This is the case e.g. in the reaction 
77 — > WW where the photons are not monochromatic but have a Compton-energy 
spectrum [20]. In the most general case the expressions for the optimal observables 
and the covariance matrix are somewhat complicated, basically because they contain 
the Jacobian determinant of the parameter transformation. However simplifications 
occur in various special cases, for example if the Jacobian determinant is a constant 
in phase space and therefore cancels in certain ratios and integrals. Using our results 
one is able to study the best statistically achievable sensitivity to a set of parameters 
in a reaction given a certain event number, if there are ambiguities in the reconstruc- 
tion of the final state. Such studies give important information on the capabilities 
of future machines and on how to choose the experimental settings like polarisations 
in an optimal way for physics studies. Apart from that optimal observables are an 
ideal tool for the analysis of experimental data. In either case one often needs explicit 
analytical expressions of the optimal observables. To this end we have collected here 
the necessary formulae for the case when the final state cannot be fully reconstructed 
from the measured variables. 
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